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Abstract 

^ ' The temperature induced phase transition is investigated in the 

, one-component scalar field (j)'^ model on a lattice by using Monte Carlo 

simulations. Using the GPGPU technology a huge amount of data is 
l/^ ■ collected that gives a possibility to determine the Linde- Weinberg low 

^ . bound on the coupling constant Aq and investigate the type of the 

phase transition for a wide interval of coupling values. It is found 
that for the values of A close to this bound a weak-first-order phase 
transition happens. It converts into a second order phase transition 
with the increase of A. A comparison with analytic calculations in con- 
^ ' tinuum field theory and lattice simulations obtained by other authors 

■ IS given. 
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1 Introduction 



The temperature induced phase transition in the 0{N) scalar field model 
with a spontaneous symmetry breaking (SSB) has a long history of inves- 
tigations. It was studied either by analytic methods of the quantum field 
theory or in lattice simulations (see Refs. [H |21 13] and references therein). It 
was recently observed by analytic calculations within the perturbation theory 
(PT) in the daisy, super daisy and some type beyond resummations ^ that 
a phase transition of the first order could occur in the 0(l)-model. How- 
ever, the lack of the expansion parameter happens near the phase transition 
temperature T ~ T^, for various kind resummations. So, it is impossible to 
draw a reliable conclusion about the transition type even for small values 
of the coupling constant A. In Ref. ^ some extended kind of resummations 
were used for the 0(iV)-model, and a phase transition of the second order 
was determined independently of the coupling value. Analogous results have 
been obtained in Monte Carlo (MC) simulations on a lattice. As a result, 
nowadays the general believe is that the phase transition is of the second 
order and the PT fails in this problem. However, in the 0(A^)-models, the 
results of PT calculations coincide with the lattice MC ones in the limit of 
N ^ oo, only HIE]. 

Recently, a new powerful computational platform - General Purpose com- 
puting on Graphics Processing Units (GPGPU) technology - has been put 
in force (Qj [10] that gives a possibility to generate extremely large amount 
of MC data. Therefore the accuracy of calculations can be essentially in- 
creased and it becomes possible to shed light upon hidden peculiarities and 
details of different processes of interest studied by MC simulations. One of 
such unsolved problems is the kind of the temperature phase transition in 
the 0(l)-model for small values of A and the reliability of the PT results. 
This is because there are no estimates what coupling values should be consid- 
ered as small ones. As a rule, the small values are chosen to be of the order 
A ~ 0.1—0.01. However, to make the correct choice some physical motivation 
is needed. One of possible reasons is the so-called Linde- Weinberg bound on 
the scalar field mass ^ |8j . Many years ago these authors observed that in 
models with the negative mass squared, < 0, the SSB did not happen for 
the coupling values below some scale, A < Aq. The actual value Aq depends 
on the mass parameter entering the Lagrangian. So, it is natural to consider 
the values of the coupling A ~ Aq as small ones. These values appear to be 
much smaller than the values mentioned above. 
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In the present paper we investigate the temperature induced phase tran- 
sition in the 0(l)-model in a wide interval of the couphng constant A using 
the GPGPU technology. We obtain the Linde- Weinberg bound Aq ~ 10~^. 
Then we compare the MC simulations with the hot and cold starts. In a nar- 
row interval just above the Linde- Weinberg bound, Aq < A < Ai ^ 10~^, an 
order parameter shows a hysteresis behavior near the phase transition tem- 
perature. This behavior means a phase transition of the first order. With 
the further increasing of A the hysteresis behavior becomes less pronounced 
and disappears at all reflecting a phase transition of the second order. 

The paper is organized as follows. In sect. 2 we describe the model and 
its realization on a lattice. In sect. 3 a necessary information on the MC 
simulations is given and the obtained results are adduced. Sect. 4 summarizes 
the results. 



2 The model 

In order to construct a self-consistent lattice version of the 0^-model we 
start with quantum field theory in continuous space. The thermodynamical 
properties of the model are described by the generating functional 

Z = JOif e-^[^], (1) 

where ip is a real scalar field, and the action is 

S = Jdx (^^dM^)dM^) - IrnMxf + ^Vix)'^ • (2) 

The standard realization of generating functional in MC simulations on 
a lattice assumes a space-time discretization and the probing random values 
of fields in order to construct the Boltzmann ensemble of field configura- 
tions. Then any macroscopic observable can be measured by averaging the 
corresponding microscopic quantity over this ensemble. However, the direct 
lattice implementation of ([I]) encounters an evident problem: the field (p is 
distributed uniformly in the infinite interval {—oo, oo) and there is no ran- 
dom number generator to simulate it. Of course, one can try to cut the 
interval off, but such a procedure requires the knowledge of characteristic 
scales of the model in order to separate the unimportant tails from the inter- 
val of physically important values. Instead, we prefer to rewrite the initial (j)^ 
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model in continuum space-time in the form allowing a further self-consistent 
lattice realization. 

First we introduce one-to-one transformation ^p{U) to a new field variable 
U{x) defined in the finite interval (0, 1). Let U = 0.5 corresponds to 99 = 
and (p{U) = — (p{l — U). The generating functional in terms of U reads 



(3) 



The Jacobian can be included in the action as det A = exp(Trlog^). Then 
the new field U can be easily realized by a uniform random number generator. 

The second step is the space-time discretization. For MC simulations 
we introduce a hypercubic lattice with hypertorous geometry. We use an 
anisotropic lattice with the spatial and temporal spacings and at — as/ C, 
C > 1. The scalar field is defined in the lattice sites. As a result, the 
generating functional becomes 



Z= J{dU{x)ex^ 

•' X 

where ^'{U) = dip/dU and 



-is[ip{U{x))]-Y.\ogip'[U{x)] 



(4) 



X '5 



^2 d^U{x)d^U{x) rr^ 
2 



■ (5) 



The lattice forward derivative is defined as usually by the finite difference 
operation 



d^U{x) ^ 



U{x + Gufl) — U {x) 



(6) 



where is the lattice spacing in the /i direction, ft is the unit vector in the 
direction indicated by /x. 

In the case of pure condensate field the action is determined by the po- 
tential 



V[U] 



Oj c f 1 2 2 fr rl ^ 4r 



(7) 
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This potential is topologically equivalent to the potential V{ip) = 
+Xip^/4:. It has one local maximum at t/ = 0.5 and two symmetric global 
minima at Uq and 1 — Uq. The spread between the values of the potential at 
the local maximum and the global minima is 

^'-'°«^-fH'"V|C/ol + ^.l£/„l). (8) 

The quantities Uq and Ay play a crucial role in MC simulations. Being 
equivalent in theory, different choices of these parameters can produce dras- 
tically different results in actual simulations. The reason is the finite number 
of simulations in an actual computer experiment. If rare but physically im- 
portant events could be missed, then the MC algorithm will not converge to 
the Boltzmann ensemble of configurations. In case of ~ 10^ iterations all the 
important probabilities have to be greater than 10~^. 

Considering the phase transition, one must guarantee that the MC algo- 
rithm meets the field values compatible with both the phases to choose. If 
Uq 0.5, then the broken phase can be missed since the corresponding field 
values are extremely rare events. On the other hand, in the limit t/o — )• (or 
Uq — )■ 1) the unbroken phase is washed out. It is also important to ensure a 
finite probability of the transition between those field values. The acceptance 
of non-zero condensate values of the field is ruled approximately by exp (Ay) 
at each lattice site. If Ay ^ 1, then the unbroken phase never occurs in 
actual simulations. If Ay — )■ 0, then there is no broken phase. To study the 
phase transition in the model, we choose the following conditions: 

Uo = 0.25, Ay = 1. (9) 

Thus, the half of generated field values will be between the global minima 
of the 'effective' potential, and no phase will be accidentally missed. The 
probability to prefer condensate or non-condensate values will be of order 
~ 0.5 ensuring the fast convergence of MC algorithm. Of course, the choice 
on]) is not optimal for temperatures far away from the critical temperature. 
To satisfy two conditions we use a convenient two-parameter function, 

^[U] = m^arctanh [r/(2f/ - 1) + (1 - r/)(2f/ - 1)=^] (10) 

with ^ > and r] > 0. The values of ^ and rj have to be found as the solution 
of equations dV /dU\u=Uo — &■ These equations can be written as 
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A 



where 2 is a dimensionless parameter of the model, 

1 - r]{2U - 1) - {1 - r]){2U - 1)^ 



JC[U] 

g[u] 



1 + ri{2U - 1) + (1 - ri){2U - 1)3 ' 
flC"[U]IC[U] 



iic'[u]y 



1 log/C[f/], 



(12) 

(13) 
(14) 
(15) 



where the primes denote derivatives with respect to U. Eq. f|TT]) gives ri{z), 
then ^ can be found from f lT2|) . There is no physical solution for z < Zmin- 
This forbidden interval corresponds to low temperatures which cannot be 
reached within the chosen parametrization. Finally the lattice action is 



Spix)] 



EE 



+ E 



■ ( ]C'[Uix)] y ( U{x + a,fx)-Uix) Y 

U^[lC[U{x)]) [ a J as ) 



-g[U{x)]~Y\og^lC[U{x)]+Y^z\og^lC[U{x)] + Vo 



Y 



21og2/C[f/o] ' 
—Ay — log(2m,^?7). 



(16) 
(17) 



The constant part of the action Vq is completely unimportant for calculations 
and can be omitted, since MC algorithm is based on the difference between 
the actions of modified and initial field configurations. 

By varying C, it is possible to change a*, while keeping fixed. Conse- 
quently the temperature T ( can be changed continuously at the fixed a^. 
At low temperatures the global discrete symmetry is expected to 

be broken due to a non-zero field condensate. The field condensate, (p, can 
be measured directly as the average of (f{x) over the lattice. As an exam- 
ple, in FigJT] we plot \ip\ in the units of the classical condensate m/y/^ for 
A = 5 ■ 10~^ and 16^ lattice. The lower values of ( and z correspond to lower 
temperatures. One can see the evident phase transition with the field con- 
densate growing with the temperature decreasing. Since the clear positive 
or negative values of ip appear in lattice configurations in the broken phase. 
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Figure 1: \ip\ in the units of classical condensate m/\/~\ for A = 5 ■ 10 ^ on 
lattice 16^ 

we conclude about the absence of domains and will use |<^| in plots in what 
follows. 

The field condensate is an obvious order parameter vanishing in the high- 
temperature phase. We can use it to determine the type of the phase tran- 
sition. In case of the first-order transition the overheated and supercooled 
states are possible. So, the MC simulations with the hot and cold starts 
have to lead to different phases near the critical temperature. Combining 
the MC simulations for the hot and cold starts we will see a hysteresis plot. 
This exfoliation in the vicinity of the critical temperature has to be observed 
independently of the direction in the (C — z) plane. 

We consider in details two slices of the two-dimensional function in FigJT] 
at fixed z = 0.35 and z = 0.5. We compute the field condensate with the 
hot and cold starts for different (. The hysteresis- type plots will mean a first 
order phase transition. 

3 Monte Carlo simulation results 

The exfoliation of the simulated data in the vicinity of the critical tempera- 
ture is a tiny effect. A large amount of simulation data must be prepared to 
observe it. In this regard, achieving the highest performance of the compu- 
tational hardware is a problem of great importance. To speed up essentially 
the simulation process we apply a GPU cluster of AMD/ATI Radeon CPUs: 
HD5870, HD5850, HD4870 and HD4850. The peak performance of the clus- 
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ter is up to 8 Tflops. The low-level AMD Intermediate Language (AMD IL) 
is used in order to obtain the maximal performance of the hardware. Some 
technical details of MC simulations on the ATI GPUs and review of AMD 
Stream SDK are given in Ref. [9] and references therein. 

The MC simulations are realized at hypercubic lattices up to 64^. Most of 
the obtained statistics come from the lattice 16^. We use the pseudo-random 
number generator RAN LUX in the MC kernel, but all the key results are 
checked with RAN MAR generator [10]. The lattice data are stored with the 
single precision. Updating the MC configurations are also performed with 
the single precision, whereas all the averaging measurements are carried out 
with the double precision to avoid the accumulation of errors. 

It should be noted that for a rough estimate of the parameter space in 
the early stages of the work we used a pseudorandom number generator 
XOR128 to speed up the calculations. The condensate values obtained by 
this generator coincide with the correspondent points produced with the 
help of RANMAR and RANLUX generators. However, the generator XOR128 
realized the two minima of the broken phase with an unequal frequency. 
Its unexpected behavior had been discovered in the first time, since good 
statistical properties of this generator were reported in the literature on real 
modeling before. 

The system is thermalized by passing 5000 MC iterations for every run. 
For measuring we use 1024 MC configurations separated by 10 bulk updates. 



As it was discussed in the previous section, we collect the data for the 
absolute value of the averaged field \(p\ representing the field condensate. The 
temperature dependence of \ip\ for the lattice 16"^ at ^ = 0.35 and z = 0.5 for 
( = [1.5; 2.5] is shown in FigJ21 The whole data set for every plot is divided 
into 15 bins. Different initial conditions are marked with different colors: the 
hot start is depicted in red (lower bins) and the cold start is represented in 
blue (upper bins). The mean values and the 95% confidence intervals are 
shown for each bin. Every bin contains 150 simulated points. 

As it is seen from FigJJl for A = 0.01 the temperature dependence of 
the field condensate is insensitive to the start configuration chosen. Both 
the cold and hot starts lead to the same behavior of the field condensate 
for various (. This means the phase transition to be of the second order, 
and this result is in agreement with the common opinion on the type of the 
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Figure 2: Temperature dependence of the absolute value of the averaged field 
\^\ for the lattice 16^ at z = 0.35 and z = 0.5 for C = [1.5; 2.4]. 



phase transition stated in Refs. [H [21 [5] . However, therein this value of the 
coupling is considered small one. 

Then, for smaller values of A the overheated configurations occur in the 
broken phase for the hot start, and the supercooled states can be found for 
the cold start. That is, the exfoliation of the simulated data in the vicinity of 
the critical temperature is observed for different start configurations. Such a 
hysteresis behavior corresponds to the phase transition of the first order. 

With further decreasing of A to the values of order Aq ~ 10~^ the behav- 
ior of hot- and cold-started simulations becomes completely separated and 
independent of the temperature. This effect is plotted in FigJS] Such type 
property means that the SSB does not happen even at zero temperature. The 
corresponding value Aq can be identified as the Linde- Weinberg low bound. 

Comparing the plots for two different values of z with the same set of other 
parameters in FigJJj one can see that the hysteresis behavior for z = 0.5 is 
similar to the case of z = 0.35. This is the common feature for all the tested 
values of A, but we adduce only one slide for brevity. The independence of 
z is important. It means that the observed exfoliation of simulated points 
is stable within our parametrization of the action. Thus, the introduced 
approach is reliable at the transition temperature. 
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Figure 3: The Linde- Weinberg regime for the lattice 16 at A = 10 , 5-10 
z = 0.35 and C = [1.5; 2.4]. 



4 Conclusion 

As it was discovered in MC simulations, the temperature phase transition in 
the 0(1) (f)^ model is strongly dependent on a coupling value A. There is the 
low bound Aq ~ 10~^ determining the range where the SSB is not realized. 
Close to this value in the interval 10~^ < A < 10"^ the phase transition is 
of the first order. For larger values of A the second order phase transition 
happens. These observations, in particular, may serve as a criterium for 
applicability of different kind resummations in perturbation theory. In fact, 
we see that the daisy and super daisy resummations give qualitatively correct 
results for small values of A. For larger values they become non-adequate to 
the second-order nature of the phase transition. In this case other more 
complicated resummation schemes should be used. 
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